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Abstract 

The classic central limit theorem and a-stable distributions play a key role in probability 
theory, and also in Boltzmann-Gibbs (BG) statistical mechanics. They both concern the paradig- 
matic case of probabilistic independence of the random variables that are being summed. A 
generalization of the BG theory, usually referred to as nonextensive statistical mechanics and 
characterized by the index q (q = 1 recovers the BG theory), introduces special (long range) 
correlations between the random variables, and recovers independence for q = 1. Recently, a 
g-central limit theorem consistent with nonextensive statistical mechanics was established [T] 
which generalizes the classic Central Limit Theorem. In the present paper we introduce and 
study symmetric (q, a)-stable distributions. The case q — 1 recovers the Levy a-stable distri- 
butions. 



1 Introduction 

In paper [T| a generalization of the classic central limit theorem applicable to nonextensive statistical 
mechanics [2j [3] (which recovers the usual, Boltzmann-Gibbs statistical mechanics as the q = 1 
particular instance), was presented; for reviews of nonextensive statistical mechanics see jUdJE]. 
We follow here along the lines of that paper. One of the important aspects of this generalization 
is that it concerns the case of random variables correlated in a special manner. On the basis of 
the q- Fourier transform F q introduced there (F\ being the classical Fourier transform), and the 
function 

we described attractors of conveniently scaled limits of sums of ^-independent random variables 
with a finite (2q — l)-variance El- This description was essentially based on the mapping 

F q :G q [2]^g z(q) [2], (1) 



1 q- independence corresponds to standard probabilistic independence if q = 1, and to specific long-range (global in 
physics terminology) correlations if q ^ 1. 

2 We required there 1 < q < 2. Denoting Q = 2q — 1, it is easy to see that this condition is equivalent to the 
finiteness of the Q-variance with 1 < Q < 3; see also [7|. 



where Q q [2] is the set of g-Gaussians (the number 2 in the notation will soon become transparent). 

In the current work, which consists of two parts, we will introduce and study a (/-analog of 
the a-stable Levy distributions. In this sense, the present paper is a conceptual continuation of 
paper [1J. For simplicity we will analyze only symmetric densities in the one-dimensional case. 
The classic theory of a-stable distributions C(a) was originated by Paul Levy and developed by 
Levy, Gnedenko, Feller and others (see, for instance, [9j [101 [HI 02] and references therein for 
details and history). The a-stable distributions found a huge number of applications in various 
practical studies [131 HH [13 HH [HI HH EQl EH E21 [23], confirming the frequent nature of these 
distributions. 

We introduce a class of random variables C q (a), which we call ((/, a)-stable distributions. 

_ 1+2 

Namely, we consider the symmetric densities f{x) with asymptotics / ~ C\x\ i+at?- 1 ), \x\ —* oo, 
where 1 < q < 2, 0<a<2, and C is a positive constant We establish that linear combinations 
and scaling limits of sequences of (/-independent random variables with (q, a)-stable distributions 
are again random variables with (q, a)-stable distributions. These facts justify that C q {a) form a 
class of stable distributions. To this end, we note that C q {a) for fixed q £ [1,2) and a G (0,2) 
coincides with the set of symmetric Levy distributions C sym (^), where 

a(2 - q) 

7 = 7(9, a) 



l + a(q- 1) 



However, the ((/, a)-stability holds for specifically correlated random variables ((/-independence 
exhibits a special correlation). Thus for each fixed q € [1, 2) we have the set of (/-independent (q, a)- 
stable distributions. If q = 1 then (/-independence becomes usual independence and 7(1,0) = a, 
implying £i(a) = £ syrn {a). The main purpose of the current paper is to classify (q, a)-stable 
distributions in terms of their densities depending on the parameters 1 < q < 2 (or equivalently 
1 < Q < 3, Q = 2q — 1) and < a < 2. We establish the mapping 

F q : G qL [2] ^ Q q [a], (2) 
where G q [ot\ is the set of functions {be q , b > 0, f3 > 0}, and 

~ L - 3 + Qa ,Q = 2q-l, 



i.e. 



1 + a 

1 + a 



q L -l 1 + «(<?- 1)' 

The particular case q = Q = 1 recovers q L = f^§, already known in the literature [3]. Denote 
Qi = {(Q,a) : 1 < Q < 3, a = 2}, Q 2 = {(Q,a) : 1 < Q < 3, < a < 2} and Q = Q x U Q 2 - 
Note that the case (Q, a) £ Q\ for (/-independent random variables with a finite Q-variance was 
studied in [T]. For (Q, a) S Q2 the Q- variance is infinite. We will focus our analysis namely on the 
latter case. Note that the case a = 2, in the framework of the present description like in that of 
the classic a-stable distributions, becomes peculiar. 

In Part II we study the attr actors of scaled sums, and expand the results of paper PQ to the 
region Q generalizing the mapping (pQ) into the form 

FUa) : W -> G Za{q) [a), 1 < q < 2, < a < 2, (3) 
3 Hereafter g{x) ~ h{x),x — > a, means that lim^ 



2 



where 

_ a-2(l-g) _ ag + 1 - g 

Q a {s)- a <md z a (s)- a + 1 _ g _ 

Note that, if a = 2, then £2(9) = g and 22(g) = (1 + g)/(3 — q), thus recovering the mapping JT]), 
and consequently, the result of p]. 

These two descriptions of (q, a)-stable distributions, based on mappings ([2]) and ([3]) respectively, 
can be unified to the scheme 

C(q,a)^G q (a) F ^ G qa j2) (4) 
I 

M 2 )> 

which gives the full picture of interrelations between the values of parameters q £ [1,2) and a £ 
(0,2). 



2 Basic operations of g-algebra 

We recall briefly the basic operations of g-algebra. Indeed, the analysis we will conduct is entirely 
based on the g-structure of nonextensive statistical mechanics (for more details see \2A\ |4] and 
references therein). To this end, we recall the well known fact that the classical Boltzmann- 
Gibbs entropy Sbg = —J2iPi^ n Pi satisfies the additivity property. Namely, if A and B are two 
independent subsystems, then Sbg(A + B) = Sbg(A) + Sbg(B). However, the g-generalization of 

l-V p q 

the classic entropy introduced in [2] and given by S q = — 1 with q £ 1Z and S\ = Sbg, does not 
possess this property if q 7^ 1. Instead, it satisfies the pseudo- additivity (or q- additivity) [2| [3"1 [22] 

S q (A + B) = S q {A) + S q {B) + (l-q) S q (A) S q (B). 

Inherited from the right hand side of this equality, the q-sum of two given real numbers, x and 
y, is defined as x (B q y = x + y + (1 — q)xy. The g-sum is commutative, associative, recovers the 
usual summing operation if q = 1 (i.e. x ©1 y = x + y), and preserves as the neutral element (i.e. 

x © 9 = x). By inversion, we can define the q-subtraction as x Q q y = l+ *^ q ^ y ■ The q-product for 

1 

x, y is defined by the binary relation x® q y = [x 1_ ' ? + y 1_ ' ? — 1]^T 9 . Here the symbol [x] + means that 

[x} + = x if x > 0, and [x]+ = if x < 0. This operation also commutative, associative, recovers the 

usual product when q = 1, and preserves 1 as the unity. The g-product is defined if x l ~ q + y 1 ~ q > 1. 

1 

Again by inversion, it can be defined the q-division: x <2> q y = [x q — y q + 1) 1-9 . 



3 (/-generalization of the exponential and cyclic functions 

Now we introduce the g-exponential and q- logarithm |24j . which play an important role in the 

nonextensive theory. These functions are denoted by e q and ln q x and respectively defined as 
1 

e q = [l + (l — q)x]]T 9 and ln q x = x 1 (x > 0). The entropy S q can be conveniently rewritten 
in the form S q = J2iPi 1% j.. 

Let us mention now the main properties of these functions, which we will use essentially in 
this paper. For the g-exponential the relations e^ qV = e q e^ and e q +y = e q © g e y q hold true. 
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These relations can be written equivalently as follows: ln q (x ® q y) = ln g x + ln g y@, and ln ? (xy) 
(Inqx) ©g (ln 9 y). The (/-exponential and (/-logarithm have the asymptotics 



e x q = l+x + |x 2 + o(x 2 ), x^O, (5) 



and 



ln g (l + x) = x-|x 2 + o(x 2 ), x^O, (6) 
respectively. The (/-product and (/-exponential can be extended to complex numbers z = x + iy (see 

In addition, for q ^ 1 the function e q can be analytically extended to the complex plain except 
the point zq = —1/(1 — q) and defined as the principal value along the cut (-co, Zq). If q < 1, then, 
for real y, \ef\ > 1 and [e*"| ~ 1^(1 + y 2 )^) , y -► oo, with K q = (1 - g)V(i-9). Similarly, if 
g > 1, then < \e q y \ < 1 and |e^| — > if \y\ —> oo. 

Lemma 3.1 Let A n (q) = Ilfc=o a k(q)i where a^iq) = q — k(l — q). Then the following power series 
expansion holds 

e * = i + z + z 2 y AlM^ z n \ z \ < 

' ^ (n + 2)! l 1 "^ 

Corollary 3.2 Lei / g = ( — 1/(1 — g|, 1/|1 — q\). For arbitrary real number x £ I q the equation 

Jx _fi_2V (-l)"^2»(g) 2ni , -r _jf (-l)M 2n+i (g) 2w+1 

holds. 

Define for x £ I q the functions (/-cos and (/-sin by formulas 

cos 9 (x) = l-x E (2n + 2) , x , (7) 

and 

^ {-lYA 2n+1 {q) 



sm q \x) = x — x 



2 \~±) fWW Jn+l /o\ 

In fact, cos g (x) and sin g (x) is defined for all real x by using appropriate power series expansions. 
Properties of (/-sin, (/-cos, and corresponding (/-hyperbolic functions, were studied in [27J. Here we 
note that the (/-analogs of the well known Euler's formulas read 

Corollary 3.3 (i) e q x = cos q (x) +i sin g (x); 
(ii) cosg(x) = q 2 q ; 
(Hi) sin q (x) 



p ix_ p -ix 



Lemma 3.4 The following equality holds: 



cos g (2x) = e 2 2q l _? )x - 2 sin^.^x). (9) 



'This property reflects the possible extensivity of S q in the presence of special correlations [30] I3T1 l32l [33] . 
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Proof. The proof follows from the definitions of cos q (x) and sin q (x), and from the fact that 
(eg) 2 = e 2 f +q y 2 (see Lemma 2.1 in [lj). ■ 

Denote ^ q (x) = cos g 2x — 1. It follows from Equation © that 

= (4-"i^ 2 " !) " 2 s <-iW- ( 10 ) 



Lemma 3.5 Lei q > 1. T/ien we /iaue 

1. -2 < < 0; 

2. ^g(x) = -2 g x 2 + o(x 3 ), x -> 0. 

Proof. It follows from (|10p that ^(x) < 0. Further, sin g (x) can be written in the form (see 
[27]) sin g (a;) = jO g (x) sin[^(a;)], where p q (x) = (e q l q ^ x ) 1//2 and tp q (x) = arcta ^(i~g) :E _ xhis yields 
^q( x ) > — 2 if g > 1. Using the asymptotic relation ([5]), we get 

4 { q~l )x2 - 1 = 2(1 - q)x 2 + o(x 3 ), x -> 0. (11) 

It follows from ([8]) that 

- 2 sin| g _ 1 (a?) = -2x 2 + o(x 3 ),x -» 0. (12) 
The relations (|1U|) . and (|12|) imply the second part of the statement. ■ 

Remark 3.6 It is not hard to verify that in the case q > 1 for x > (q — 1) _1 the representation 
q M ) l { (l-g) 2 x (l-g) 4 x 2 ^ (n + 2)!(g-l) 2 ™ V ) 



follows from Lemma \3.1l 

4 g-Fourier transform for symmetric functions 

The q- Fourier transform for q > 1, based on the g-product, was introduced in [1] and played 
a central role in establishing the (/-analog of the standard central limit theorem. Formally the 
q- Fourier transform for a given function f(x) is defined by the formula 

/oo 
eft® q f(x)dx. (13) 

For discrete functions fk,k = 0, ±1, this definition takes the form 

oo 

F q [f}(0= E e f ( 14 ) 

k=— oo 

In the future we use the same notation in both cases. We also call (|13p or (|14p the q- characteristic 
function of a given random variable X with an associated density f(x), using the notations F q (X) 
or F g (f) equivalently. 
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It should be noted that, if in the formal definition (fT3|) . / is compactly supported, then integra- 
tion has to be taken over this support, although, in contrast with the usual analysis, the function 
e q x ^ ® q f{x) under the integral does not vanish outside the support of /. This is an effect of the 
g-product. 

The following lemma establishes the relation of the g-Fourier transform without using the q- 
product. 

Lemma 4.1 The q-Fourier transform can be written in the form 

/( x )e^(/(»))'- dx. (15) 

-oo 

Remark 4.2 Note that, if the q-Fourier transform of a given function f(x) defined by the formal 
definition in \13\) exists, then it coincides with the expression in U5\) . The q-Fourier transform 
determined by the formula [T5\) has an advantage when compared to the formal definition: it does 
not use the q-product, which is, as we noticed above, restrictive in use. 

Further to the properties of the q- Fourier transform established in PQ, we note that, for sym- 
metric densities, the assertion analogous to Lemma 14.11 is true with the q-cos. 

Lemma 4.3 Let f(x) be an even function. Then its q-Fourier transform can be written in the 
form 

/oo 
f(x)cos q (x^[f(x)}i- 1 )dx. (16) 
-oo 

Proof. Notice that, because of the symmetry of /, 

/OO /"OO 
ef e ®« f(x)dx = / e~ ix ^ (S> q f(x)dx. 
-OO J —oo 

Taking this into account, we have 

Fq[m = 2j 00 ( e " ? ® q f{x) + ^ ® 9 dX ■ 

Applying Lemma 14.11 we obtain 

F q [f](0 = / m~ V s dx, 

J —oo & 

which coincides with (|16p . ■ 

Further, denote FL q ^ a = {/ G L\ : fix) ~ C\x\ i+<*(?-i) , \x\ — > oo}. It is readily seen that 
<M<7, a) = i+afo-i) > 1 for a11 a G (°> 2 ) and Q G I 1 ) 2 )- Moreover, c/)(q,a)(2q - 1) < 3 for all 
a G (0,2) and q G [1,2), which implies cr|„ x (/) = oo. 

Lemma 4.4 Let f(x), x G R, be a symmetric probability density function of a given random vari- 
able. Further, let either 

(i) the [2q — l)-variance 0"|„_i(/) < ooH (associated with (2q — 1, a) G Qi)> or 

(ii) f{x) G H Qi a, where {2q - l,a) G Q 2 . 

5 In [I] we did not require the condition for f(x) to be symmetric if tJ2q-i(f) < oo. 
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Then, for the q-Fourier transform of f(x), the following asymptotic relation holds true: 

F q [f](C) = l-V q M a + om a ),^0, (17) 

where 



( q 2 



f^q,a — < 



if a = 2 ; 



(18) 



2 2 -"(l+g(g~l))C roo - ggfa) 



2-9 



Jo°°^#%> if(2q-l,a)€Q 2 



w l thv 2q ^(f) = JZ [f(x)] 2q - 1 dx . 



Remark 4.5 Stable distributions require to be positive. We have seen (Lemma \3. 5\) that if 
q > 1, i/ien ^(x) < (Vioi being identically zero), which yields Hq t a 

> 0. 



Proof. First, we assume that a = 2. Using Lemma 14.11 we have 

/oo roo 
{^)®qf{x)dx= / /(x)cOS g K[/(x)]^ 1 )x. (19) 
-oo J —oo 

Making use of the asymptotic expansion ([5]) we can rewrite the right hand side of (|19p in the form 

/OO , . 

/(x) (l + ixt[f{x)]i- 1 - q/2x 2 eif(x)] 2 ^) dx + (£ 3 ) = 

1 - (g/2)? 2 <r 2 V 1 ^-i + o(£ 3 ), £ - 0, 

from which the first part of Lemma follows. 

Now, assume (2q — 1, a) £ Q2- We apply Lemma 14731 to obtain 

/oo 
/(x)[cos 9 (xe[/(x)]^ 1 )-l]dx = 
-oo 

2/ f{x)*q{ o )dx + 2 f{x)tyq{ )dx, 

Jo 2 Jn 2 

where N is a sufficiently large finite number. In the first integral we use the asymptotic relation 
\l/(§ ) = — \x 2 + o(x 3 ), which follows from Lemma and get 

2^/(x)* g (^^)dx = 

/•JV 

- / x 2 fi-\x)dx + 0(f), C - 0. (20) 
Jo 

In the second integral taking into account the hypothesis of the lemma with respect to /(x), we 
have 

•/AT ^ ./AT X l + a(l-q) ^ 

We use the substitution 

2 ~i 2y 

x l + a(q-l) = _2. 
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in the last integral, and obtain 

f x£[f(x))i- 1 



JN 2 



q Jo y 

Hence, the obtained asymptotic relations (f20|) and (f2"T|) complete the proof. 



5 Weak convergence of correlated random variables 

This section we start with introduction of the notion of q-independence in particular case. More gen- 
eral definition and some examples are given in pQ. We will also introduce two types of convergence, 
namely, q-convergence and weak q-convergence and establish their equivalence. 

Definition 5.1 Two random variables X and Y are said to be q-independent® if 

F q [X + ¥](£) = F q [X] (£) ® 9 F g [Y](£) . (22) 

In terms of densities, the relation (|22p can be rewritten as follows. Let fx and fy be densities of 
X and Y respectively, and let fx+Y be the density of X + Y. Then 

/oo 
ej* ® g /x+y (x)dx = F q [f x m ® g ^[/y](6- (23) 

For o = l the conditions (|22|) turns into the well known relation 

^[/x * /y] = F[f x ] ■ F[f Y ] 

between the convolution (noted *) of two densities and the multiplication of their (classical) Fourier 
images, and holds for independent X and Y. If q ^ 1, then g-independence describes a specific 
correlation. 

Definition 5.2 Let Xi,X2, ...,X^, ... be a sequence of identically distributed random variables. 
Denote Y^ = X\ + ... + X^. By definition, the sequence Xn,N = 1,2, ... is said to be q-independent 
(or q-i.i.d.) if for all N = 2,3, the relations 

F g [Y N ](0 = F q [Xx](0 ® q ... ® q F q [X N ](0 (24) 

hold. 

Definition 5.3 A sequence of random variables Xn, N = 1,2, is said to be q-convergent to a 
random variable X^ if limjv^oo F q [XN](C) = i ? g [X oc ](^) locally uniformly in £. 

Evidently, this definition is equivalent to the weak convergence (denoted by "=>") of random 
variables if q = 1. For q ^ 1 denote by W q the set of continuous functions (p satisfying the condition 

100*01 < C(l + \x\y^, xeR. 

6 We assume X and Y to have the zero g-means. For the definition of g-independence for random variables with 
non-zero g-means see [I] 
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Definition 5.4 A sequence of random variables X^j with the density /jy is called weakly q-convergent 
to a random variable with the density f if f R f]y(x)dm q — ► J Rd f(x)dm q for arbitrary measure 
m q defined as dm q (x) = 4> q (x)dx, where (j) q G W q . We denote the q-convergence by the symbol 4>. 

Lemma 5.5 Let q > 1. Then Xjy => Xq yields Xjy => Xq. 

The proof of this lemma immediately follows from the obvious fact that W q is a subset of the 
set of bounded continuous functions. 

Recall that a sequence of probability measures is called tight if, for an arbitrary e > 0, there 
is a compact K € and an integer N* such that [i]\[(R d \K e ) < e for all N > N*. 

Lemma 5.6 Let 1 < q < 2. Assume a sequence of random variables Xn, defined on a probability 
space with a probability measure P, and associated densities /at, is q-convergent to a random variable 
X with an associated density f. Then the sequence of associated probability measures [In = PiX^ 1 ) 
is tight. 

Proof. Assume that 1 < q < 2 and Xn is a (/-convergent sequence of random variables with 
associated densities /at and associated probability measures ^n- We have 

i J\l - F g [f N }(0)dti = i jf* (1 - J R f N e l q U - 1 dx)di = 

I J R{1 " Ul£ ) dm{x) - (25) 



It is not hard to verify that 



! fR 2sm^(Rx(2-q)t) 

/ e «#df = 2 -±- . (26) 

RJ-R q Rx(2-q)t v ' 



It follows from ([25]) and (f26l) that 



i r R roo ( sin^(x(2 - q)Rf q N X ) 

R I J 1 - F < [f » md t ='L (' - "r,(2-^ 1 ^ (2?) 

Since 1 < q < 2 by assumption, > 1 as well. It is known [271 US EH] that for any q > 1 the 
properties sin q i(x) < 1 and (sin q i (x))/x — ► 1, x — > hold. Moreover, (sin / (x))/x < 1, Vx G i?. 

Suppose, lim\ x \^ >00 \x\f'^ 1 = Lat, iV > 1. Divide the set {N > iVo} into two subsets A = {Nj > 
Nq : Lj\r. > 1} and B = {N^ > Nq : Ln < 1}- If N G A, since sin^_ < 1, there is a number a > 
such that 

1 



— d/j, N (x) 



- f (1 - F v [f N ](€))d£ > 2 f (1 

RJ-R K ? ^ UWM - 4|>aV ^1(2-^ ; 

> CfjL N (\x\ > a), C > VJVei, 

for i? small enough. Now taking into account the g-convergence of Xjy to X and, if necessary, 
taking R smaller, for any e > 0, we obtain 



fi N [\x\ 



>a)<-^j\l-F q [f ](0)d^<e, ViV G A. 
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If N E B then there exist constants b > 0, 5 > 0, such that 

fN{x) < < i±A | X | > 6 ,viV E B. 

\x\i- 1 \x\i~ 1 

Hence, we have 

m (\ x \ > b) = J f N {x)dx < (1 + 6) J ——j-, N e B. 

J\x\>b J\x\>b | | <? — i 

Since, l/(q — l) > 1, for any e > we can select a number b e >b such that //jv(|x| > b e ) < e, N E -B. 
As far as U -B = {iV > Ao} the proof of the statement is complete. ■ 
Further, we introduce the function 

D q (t) = D g (t; a) = tef' 1 = t(l + i(l - q)at q - 1 y^ , (28) 

defined on [0,1], where 1 < q < 2 and a is a fixed real number. Obviously D q (t) is continuous 
on [0, 1] and differentiable in the interval (0, 1). In accordance with the classical Lagrange average 
theorem for any ti,t2, < t\ < t 2 < 1 there exists a number t#, t\ < t* < t 2 such that 

D q (h) - D q (t 2 ) = D' q (U){h - t 2 ), (29) 

where D q means the derivative of D q (t) with respect to t. 

Consider the following Cauchy problems for the Bernoulli equation 

y'-\y = ^ 1 ^y«,y(o) = o, (30) 

It is not hard to verify that D q {t) is a solution to the problems (|30p . 
Lemma 5.7 For DJt) the estimate 

\D q (t;a)\ < C(l + \a\y^, t E (0,1], a E R 1 , (31) 
holds, where the constant C does not depend on t. 
Proof. It follows from fl2SJ and ([30]) that 

\y\t)\ < t' l \y + ia(l - q)y«\ = \ef ^ + ia(l - g)(ef ">| = 
|1 + io(l - g)^ 1 ! - ^ < (7(1 + |a|)~«^, t E (0, 1]. 



Theorem 5.8 Let 1 < q < 2 and a sequence of random vectors Xn be weakly q-convergent to a 
random vector X . Then Xjy is q-convergent to X. 

Proof. Assume Xn, with associated densities fx, is weakly g-convergent to a X, with an 
associated density /. The difference ^ [/tv] (<!; ) — •T'gLfivKO can be written in the form 



^[/tf](0"-Wiv](0= / (D q {f N (x))-D q {f(x)))dx, (32) 

JR d 

where D q (t) = D q {t; a) is defined in ([28]) with a = x£. It follows from ([29]) and §B) that 

9 



WjvKO " -WivKOI < C / |(1 + |x|)-— (/jv(x) - /(x)) |dx, 
which yields ^[/jv](0 - ^,[/iv](0 for all £ E ■ 
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Theorem 5.9 Let 1 < q < 2 and a sequence of random vectors X/v with the associated densities 
/at is q-convergent to a random vector X with the associated density f and .Fg [/](£) is continuous 
at £ = 0. Then weakly q-converges to X. 

Proof. Now assume that /at converges to / in the sense of ^-convergence. It follows from 
Lemma [5.61 that the corresponding sequence of induced probability measures [In = -PpsT^ 1 ) is 
tight. This yields relatively weak compactness of the sequence hn- Theorem 15.81 implies that each 
weakly convergent subsequence {/J-Nj} of /i^ converges to n = P{X~ l ). Hence, /Uat =4> ^, or the 
same, Xn => X. Now applying Lemma 15.51 we complete the proof. ■ 



6 Symmetric (g, a)-stable distributions. First representation 

In this section we introduce symmetric (q, a)-stable distributions and give the description based 
on the mapping ([2]). In accordance with this description q takes any value in [1,2), however we 
distinguish the cases a = 2 and a G (0, 2). 



Definition 6.1 A random variable X is said to have a (q,a)-stable distribution if its q-Fourier 
transform is represe 
stable distributions. 



transform is represented in the form e q , with (3 > 0. We denote by C q (a) the set of all (q,a) 



Denote Q q (a) = {be q ma , b > 0, > 0}. In other words X G C q {a) if F q [f) G Q q {a) with 6=1. 
Note that if a = 2, then Q q {2) represents the set of g-Gaussians and C q {2) - the set of random 
variables whose densities are g*-Gaussians, where q* = (3q — 1)/(1 + q). 

Proposition 6.2 Let q-independent random variables Xj G C q (a),j = 1, ..,m. Then for any con- 
stants oi, a m , 

m 

^ajXj G C q (a). 
i=i 

Proof. Let 

F q [X j }(0 = e-^\Z\<*,j = l,...,m. 

Using the properties e x q ® q e y q = e q +y and F q [aX](£) = F q [X\{a?~ q ^), it follows from the definition 
of the ^-independence that 

m m 

F q [J2a j X j ] = e-m^p = j :Pjln l^- q) 

3=1 3=1 



Remark 6.3 Proposition 1 6. S\ justifies the stability of distributions in C q {a). Recall that if q = 1 
then q-independent random variables are independent in the usual sense. Thus, if q = 1, < a < 2, 
then £\(a) coincides with symmetric a-stable Levy distributions C sym {a). 

Further we show that the g-weak limits of sums 

1 



J N 



(X 1 + ... + X N ),N = l,2,.. 



s N (q,a) 

as iV — > oo, where S]y(q, «), N = 1, 2, are some reals (scaling parameter), also belong to C q {ot). 
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Definition 6.4 A sequence of random variables Zjy is said to be q-convergent to a (q, a) -stable 
distribution, if lim^ F q [Z € Gq{ot) locally uniformly by £. 

Theorem 1. Assume (2q — l,a) £ Qi- Let Xi, X2, Xn , ... be symmetric q-independent 

random variables and all having the same probability density function f(x) € H qt0l . Then Z^, with 
1 

sj\f(q,a) = (fig^N) a ( 2 -9) , is q-convergent to a (q, a) -stable distribution, as N — > 00. 

Remark 6.5 By definition Q2 excludes the value a = 2. The case a = 2, in accordance with the 
first part of Lemma \4.4\ coincides with Theorem 2 of [If. Note in this case C q {2) = Q q *{2), where 

* _ 3g-l 
q ~ <7+l • 

Proof. Assume (Q,a) G Q%. Let / be the density associated with X\. First we evaluate 
F q {X\) = F q (f(x)). Using Lemma I4T41 we have 

F q [m) = 1 - M9ia ier + o(\z\ a u - 0. (33) 

Denote Yj = N~~ Xj,j = 1,2,.... Then Zjy = Y\ + ... + Yjv- Further, it is readily seen that, 
for a given random variable X and real a > 0, the equality F q [aX](l;) = F q [X](a 2 ~ q £) holds. It 
follows from this relation that F q (Yj) = F q [f](-^ — ), j = 1, 2, ... Moreover, it follows from the 
^-independence of Xi,X2, and the associativity of the g-product that 

F q [Z N ](0 = F q [f]((^ a N)-^0^-^ q F q [f}(^ q ,aN)-^0 (N factors). (34) 

Hence, making use of the expansion ([6]) for the (/-logarithm, Eq. (134ft implies 

ln q F q [Z N }(0 = Nln q F q [f]((^ a N)-^0 = Nln q (l - ^ + G (^)) = 

- |£| a + o(l), N -> 00, (35) 

locally uniformly by £. 

Hence, locally uniformly by £, 

lim F g (Z 7V ) = e-l«l Q eg q (a). (36) 

iV— >oo 

Thus, is g-convergent to a (q, a)-stable distribution, as iV — > 00. ■ 

This theorem links the classic Levy distributions with their q^-Gaussian counterparts. Indeed, 
in accordance with this theorem, a function /, for which 

/ ~ c , /x( Q+1 )/( 1+Q ^ 1 », \x\ -► 00, 

is in C q (a), i.e. F q [f](^) £ <5 g (a). It is not hard to verify that there exists a g^-Gaussian, which 
is asymptotically equivalent to /. Let us now find q^. Any (/^-Gaussian behaves asymptotically 
C\/\x\ r ' = C2/\x\ 2 /( q <* Cj = const, j = 1,2, i.e. 77 = 2/(q^ — 1). Hence, we obtain the relation 

a + 1 2 



1 + a(q - 1) g£ 
Solving this equation with respect to we have 



(37) 



q L a = 3 -^,Q = 2q-l, (38) 

linking three parameters: a, the parameter of the a-stable Levy distributions, q, the parameter 
of correlation, and q^, the parameter of attractors in terms of g^-Gaussians (see Fig. [2] (left)). 
Equation (|38p identifies all (Q, a)-stable distributions with the same index of attractor G„l (See 
Fig. H]), proving the following proposition. 
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Proposition 6.6 Let 1 < Q < 3, (Q = 2q — 1) and < a < 2. TTien all distributions X 6 C q {a), 
where the pairs (Q, a) satisfy the equation 

3 + Qa t 
a + 1 

/iai>e i/ie same attractor asymptotically equivalent to q^- Gaussian. 




Figure 1: All pairs of (Q,a) on the indicated curves are associated with the same g^-Gaussian. Two curves corresponding 
to two different values of q^ do not intersect. In this sense these curves represent the constant levels of q^ or r\ = 2/{q^ — 1). 
The line r\ = 1 joins the points (Q, a) = (1, 0.0 — 0) and (3 — 0, 2); the line r) = 2 joins the Cauchy distribution (noted C) with 
itself at (Q, a) = (1, 1) and at (2, 2); the r\ = 3 line joins the points (Q, a) = (1, 2.0 — 0) and (5/3, 2) (by e we simply mean to 
give an indication, and not that both infinitesimals coincide). The entire line at Q = 1 and < a < 2 is mapped into the line 
at a = 2 and 5/3 < q* < 3. 

In the particular case Q = 1, we recover the known connection between the classical Levy 
distributions (q = Q = 1) and corresponding g^-Gaussians. Put Q = 1 in Eq. (|38|) to obtain 

Qa = JTT^> < a < 2. (39) 

When a increases between and 2 (i.e. < a < 2), q^ decreases between 3 and 5/3 (i.e. 5/3 < 
q% < 3): See Figs. [2] (left) andEJleft). 

It is useful to find the relationship between rj = L _ , which corresponds to the asymptotic 
behaviour of the attractor depending on (a,Q). Using formula f)38h . we obtain (Fig. right) 



2(a + l) 

" = 2 + a( Q-iy <40) 

Proposition 6.7 Zei X E £g(a), 1<Q<3, 0<a<2. TTien the associated density function fx 
has asymptotics fx{x) ~ \x\ v , \x\ — > oo, where n = n(Q,a) is defined in (Jty . 

Remark 6.8 If Q = 1 (classic Levy distributions), then n = a + l, as is well known. 
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Figure 2: q^ = q* and r\ as functions of variables (Q,a). 

Analogous relationships can be obtained for other values of Q. We call, for convenience, a 
(Q, a)-stable distribution a Q-Cauchy distribution, if its parameter a = 1. We obtain the classic 
Cauchy-Poisson distribution if Q = 1. The corresponding line can be obtained cutting the surface 
in Fig. [2] (right) along the line a = 1. For Q-Cauchy distributions we have 

9l (Q) = ^ and r? = (41) 

respectively (see Figs. [2]). 

The relationship between a and for typical fixed values of Q are given in Fig. [3] (left). In 
this figure we can also see, that a = 1 (Cauchy) corresponds to q\ = 2 (in the Q = 1 curve). In 
Fig. [3] (right) the relationships between Q (Q = 2q — 1) and are represented for typical fixed 
values of a. 




Figure 3: Constant Q and constant a sections of Fig. [2fleft). 
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